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The aim of this paper is to derive a new equation (the acoustic- drift equation (ADE)) 
describing the generation of a flow by an acoustic wave. We consider acoustic waves of 
perfect barotropic gas as the zero-order solution and derive the equation for the averaged 
flow of the first order. The used small parameter of our asymptotic study is dimcnsionless 
inverse frequency, and the leading term for a velocity field is chosen to be a purely 
oscillating acoustic field. The employed mathematical approach combines the two-timing 
method and the notion of a distinguished limit. The properties of commutators are used 
to simplify calculations. The derived averaged equation is similar to the original vorticity 
equation, where the Reynolds stresses has been transformed to an additional advection 
with the drift velocity. Hence ADE can be seen as a compressible version of the Craik- 
Lcibovich equation. In particular, ADE shows that if the averaged vorticity is absent in 
the initial data, then it will always be equal to zero. This property confirms that the 
acoustic streaming is of viscous nature. At the same time ADE can be useful for the 
description of secondary motions which take place for 'strong' acoustic streaming flow. 



1. Introduction 

Acoustic streaming represents a classical topic of fluid dynamics, see e.g. Riley (2001) 
Nyborg (1998)||Lighthill (1978a)[|Lighthill (1978)[ It usually describes the averaged flow 
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generated by a given acoustic wave. The related asymptotic analysis is usually based 
on the Navier-Stokes equations and goes up to the second-order (in amplitude) terms, 
while the terms of the third and higher orders appear as physically irrelevant. In contrary 
to this classical approach we build a rigorous asymptotic procedure based on the Euler 
equations where the first-order averaged equations (we call it an acoustic-drift equation 
(ADE)) appears from the treatment of the third-order approximation (when the primary 
acoustic wave is treated as the zero order). 

Mathematically, we suggest that velocity, density, and pressure oscillate with fre- 
quency higher than inverse characteristic times of all other co-existing motions. It can 
be expressed as a 3> 1, where a is dimensionless frequency of oscillations. To de- 
rive the averaged equation we employ the two-timing method in the form used by 



Vladimirov (2012) The resulted averaged vorticity equation is similar to the original 
equation, but surprisingly a drift velocity (instead of Reynolds stresses) plays a part of an 
additional advection velocity. The derived equation for the averaged first-order vorticity 
(ADE) formally coincides with the Craik-Lcibovich equations for an incompressible fluid 
( |Craik fc Leibovich (1976)||Leibovich (1983)||Craik (1985)HVladimirov (2012)] ), while the 



difference in our case appears as the prescribed divergence of velocity. An additional in- 
centive for our research is to make all the calculations and derivations elementary and 
free of any physical and mathematical assumptions (except the most common ones, such 
as the existence of diffcrentiable solutions) . We use only Eulerian description and Eulc- 
rian averaging operation, since they are the most transparent and allow us to avoid any 
unnecessary difficulties. 
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2. General setting 

2.1. Formulation of problem 

The governing equation for an inviscid compressible barotropic fluid (gas) with velocity 
field u* , density p* , vorticity ui* = curl *it* and pressure p* is taken in the vorticity form 

du*/dt* + [u)*,u*]* +u*dbru* =0, in V* (2.1) 
dp*/dt* + {u* ■ V)p* + p*divtt* =0 = 

where asterisks mark dimensional variables and operations, i*-timc, 
cartesian coordinates, and [•,■]* stands for the dimensional commutator (|2.11[) . The term 
with pressure has been excluded since the fluid is barotropic p* — p*(p*). In this paper 
we deal with the transformations of equations, hence the form of flow domain T>* and 
particular boundary conditions can be specified at the later stages. The related equation 
for the ratio A* = lj* j p* is 

dX*/dt* + [X*,u*]* =0, in V* (2.2) 

We accept that the considered class of (unknown) oscillatory solutions u*,p* possesses 
characteristic scales of velocity U, length L, density pjj and high frequency a* 

U, L, p*, a* > 1/T; T = L/U (2.3) 

where p$ is reference density (at the state of rest), T is a dependent time-scale. The 
dimensionless variables and frequency are 

x = x*/L, t = t/T, u = u*/U, p = p*/p* a , o- = cr*T>l (2.4) 
2.2. Notations 

The variables x — [xx, x%, 2:3), t, s, and r in the text below serve as dimensionless carte- 
sian coordinates, physical time, slow time, and fast time. The used definitions, notations, 
and properties are: 
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(i) A dimensionless function / = f(x,s,r) belongs to the class O(l) if / = O(l) and all 
partial a:-, s-, and r-derivatives of / (required for our consideration) are also 0(1). In 
this paper all small parameters appear as explicit multipliers, while all functions always 
belong to ©(l)-class. 

(ii) The class H of hat-functions (or oscillating functions with the non-zero mean) is 
defined as 

/GH: /(aj ) «,r) = /(x J * J T + 27r) (2.5) 

where the r-depcndcncc is always 27r-periodic; the dependencies on x and s are not 
specified. 

(iii) The subscripts t, t, and s denote the related partial derivatives. 

(iv) For an arbitrary / £ H the averaging operation is 



1 



T0 + 27T 



(/) = —/ f(x,s,r)dT, Vr (2.6) 

where during the r-integration we keep s = const and (/) does not depend on to- 

(v) The class T of tilde-functions (or purely oscillating functions, or fluctuations) is such 
that 

/ST: J(x,s,t)=J{x,s,t + 2<k), with (/)=0, (2.7) 

they represent a special case of hat-functions with the zero average. The class B of bar- 
functions (or mean-functions) is defined as 

/£»: / r =0, f(x,s) = (J(x,s)} (2.8) 

(vi) We introduce the tilde-integral (or the fluctuating part of an integral of a fluctuating 
function) f T as 



r = J /(», s, V)dv-^J ( f f(x, s, n) d V ) dp (2.9) 



»2tt /./j 
'0 Z; " ^0 ^^0 

The tilde-integration is inverse to T-diffcrcntiation (f T ) T = {f T ) T = f', the proof is 
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omitted. The r-derivative of a tilde-function always represents a tilde-function. However 

the r-integration of a tilde-function can produce a hat-function. {For example: let us take 

4> = 0q sin t where (f> is an arbitrary bar-function: one can see that (4>) = 0, however 

(/J" (j>(x, s, p)dp) = 4> ^ 0, unless (j> = 0.} Formula (|2.9[) keeps the result of integration 

inside the T-class. 

(vii) The unique solution of a PDE inside the tilde-class 

f T = dJ/dT = o => 7=o (2.io) 

follows from (|2.9[) . One can also see that f T = f T , (f T ) = (f T ) = 0. 

(viii) The commutator of two vector-fields / and g is 

[f,g] = (g-V)f-(f-V)g (2.11) 

The commutator is antisymmetric and satisfies Jacobi's identity for any vector-fields /, 
g, h: 

[f,9] = -[g,f], [fA9,h}] + [h,[f,g]} + [g,[h,f]} = (2.12) 

As the average operation (|2.6[) is proportional to the integration over r, the integration 
by parts yields 

([f,g T ]) = -{[f T ,g]) = -<[/ T ,5]), ff,g T }) = -{[7,9]) = -{[1,9]) (2.13) 
For any tilde- function / and bar- function g (|2.12p . (|2.13[) give 

{[l,[gT]]) = [9,V] where V = {[f 7 7})/2 (2.14) 
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3. Asymptotic procedure 

3.1. Two-timing setting and distinguished limit 

We assume that a flow has its own intrinsic slow-time scale T s \ OV! (which can be different 
from T) and consider solutions of (|2.1[) in the form of hat-functions (|2.5[) 

u* = Uu(x, s, t), p* = pop{x, s, t); with t = at, s = fl t, SI = T/T s \ ow (3.1) 

Then the use of the chain rule and the transformation to dimensionless variables give 

[Q, u] H — a;divu, = (3-2) 
a 

(u ■ V)p + -pdiv u = = 
a 

In order to keep variable s 'slow' in comparison with r we have to accept that £l/cr <C 1. 
Then eqn. p.2j) contains two independent small parameters: 

11 1 CI 

£ = Tf~z = -, £ i = Tf -~z = — ( 3 - 3 ) 

la a 1 s iowCT a 

Here we must make an auxiliary (but technically essential) assumption: after the use of 

the chain rule (|3.2j) variables s and r are considered to be mutually independent: 

T, s — independent variables (3-4) 



d 


n d > 


|9- 


1 


!h ' 


a ds j 


V - 
a 


d 


fl d > 




■i( 




a ds > 







Following the same distinguished limit argument as in (Vladimirov (2012) I we take 



Vl(a) = l/a: r = at, s = t/a (3.5) 

In this paper we will show that the first three successive approximations (with s (|3.5p ) do 
produce a valid asymptotic solution. The proof of other properties of the distinguished 
limit are similar to the one considered in Vladimirov (2010) Vladimirov (2011)] this 



proof is beyond the style and scope of this paper and is omitted. 
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Hence the governing equations are 

£} T + e[cD, u] + eujdivu + e 2 Q s = (3.6) 

p T + e(u ■ V)p + epdivu + e 2 p s = (3.7) 

A r +e[X,u] + e 2 X s = (3.8) 



where e = 1/a — >• 0. Let us look for the solutions of (|3.6[) - (|3.8p in the form of regular 
series 

oo 

{Q,u,p,\) ^Y. £k ( Q k,u k ,p k ,X k ); u t ,2 fc ,p t ,A t einO(l), k = 0, 1,2, . . .(3.9) 

fc=0 

We study the class of flows which represents acoustic waves (and the related streaming 
flows) at the state of rest: 

u = 0, p = l (3.10) 

which implies that 

uJq = 0, A = (3.11) 

The considered asymptotic series for velocity and density are 

u = u + eu\ + e 2 u 2 + ... (3.12) 
p = 1 + epi + e 2 p 2 + ... 

where tto & n d pi satisfies the classical linear acoustic equations that easily follow after 
the substitution of (|3 . 1 2|) into the Euler's equations, which are related to (|3.6[) , (|3.7p but 
not presented here for brevity. An important fact for Euler's equations is: the application 
of curl^ 1 to (|3.6[) should produce a gradient of a function (enthalpy) which must be 
0(1) if we wish to build a valid asymptotic procedure; this order of enthalpy requires 
an appropriate scale (in terms of e) for pressure. In fact, all the derivation of ADE can 
be carried ahead in terms of velocity field for the Euler equations, however in this paper 
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we exploit the simplification of calculations for the vorticity equation. It means that the 
fields 

u = V(f>o(x,t), p 1 =pi(x,t) (3.13) 

everywhere below are considered as a given solution of the classical inviscid linear acoustic 
equations. 

3.2. Successive approximations 

The substitution of (|3.9I) - (|3.13I) into (|3.6[) - (|3.8[) produces the equations of successive ap- 
proximations. Let us first consider the equation for A (|3.8|) : 

The zero-order equation is Xq t = 0. Its unique solution (|2.10|) is Ao = and u>q = 0. 
Taking into account (|3 . 11 1) we can write 

A = 0, Q = (3.14) 

which means that the zero-order flow is potential and purely oscillating. Eqn (|3.14[) leads 
to the similar first-order approximation of Ai T = which have the unique solution 

Xi = 0„ A 1= [7] (3.15) 

where a mean function remains undetermined. The second-order equations that take into 
account (|3.11j) . (|3.14j) . (|3.15[) are A 2t + [Ai, «o] = which after tilde-integration (|2.9j) yield 

A 2 = K,Ai], A 2 =0 (3.16) 

The third-order equation that takes into account (|3.11[) . (|3.14l) . (13.15[) is 

A 3r + Ala + [A 2 , So] + [Al, Ui] = 

The bar-part (|2.6[) (or the averaged part) of this equation is 

A 1s + [A 1 ,m 1 ] + ([A 2 ,S ]) =0 (3.17) 
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which can be transformed with the use of (|3.16j) and (|2.14|) into the form 

X ls + [X 1 ,u 1 + Vo}=0 (3.18) 

where the drift velocity is 

V = ([u ,u T })/2 (3.19) 

where Uq is given by p,13[) . Taking into account that Ai = Ui/po p.!9[) can be rewritten 

as 

5?ia + [5?i J 5i + V'o] = (3.20) 

which is the same as 

ZJ ls + V x (tJi x (ui +V )) = widiv (Hi + V ) (3.21) 

The next step of our consideration is to repeat a similar (but analytically simpler) asymp- 
totic procedure for the density equation p. 71) . Omitting the calculations we state that it 
yields 

div(tt 1 + y o )=0 (3.22) 

and the equation (|3.2ip is reduced to 

ZJ ls + V x (w x x («i + F )) = (3.23) 

Taking curl -1 of this equation yields 

Ux s +uJi x Hi +ZUi x V — -VIIi (3.24) 

where 111 is modified pressure. The component version of this equation (where we drop the 
subscripts '1' of the first approximation and '0' for a drift, use the summation convention, 
and introduce different modified pressure p) is: 
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New equations for a streaming flow (|3.22|1 - ([3.25|1 (which we call ADE) represent the main 
result of this paper. 



4. Discussion 

1. One can see that ADE (|3.25p represents an Euler's equation 'modified' by additional 
'advection term' and 'linear friction term' 

dV k _ 
— 5 — Uk 

OXi 

This 'effective friction' is linear in velocity u and can be positive or negative. This 
term is important for the zero-order solutions (|3.13[) with the amplitude of acoustic 
wave changing in space, similar to described in Lighthill (1978a){ Nyborg (1998) Indeed, 



with the decreasing of the amplitude Uo| of an acoustic wave (|3.13[) with the distance 
(say, with the increasing the distance from the source of sound), the value \Vq\ (|3.19[) 
is also decreasing. Therefore the 'friction term' become positive (which corresponds to 
'anti- friction' or 'negative friction'); from the physical viewpoint it should generate the 



acoustic streaming. This argument formally repeats the treatment of Lighthill (1978a) 



Nyborg (1998) However returning to the averaged vorticity equation (|3.23[) , one can see 
that if the average vorticity is zero at the initial data, then it will always be equal to 
zero. It means that the presence of a variable in space amplitude is not sufficient for 
the generation of acoustic streaming. This conclusion confirms the classical idea about 
essentially viscous nature of acoustic streaming. At the same time ADE contradicts to the 
estimations of acoustic streaming obtained without taking viscosity into consideration, 



see 



Nyborg (1998) p. 214, where a prescribed attenuated plane wave is introduced 'by 
hands' into an inviscid setting for forces. 

2. The above explanation indicates that ADE does not describe the acoustic streaming 
when the streaming velocity is small. However ADE can be relevant to 'intense' regimes 
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of acoustic streaming, including the generation of secondary vortices, similar to Langmuir 



circulations Craik fc Leibovich (1976)] Vladimirov (2012) 



3. Our approach (based on the two-timing method, distinguished limit, and the unusual 
choice of small parameter) is different from the mainstream classical theory of acoustic 
streaming. It suggests that ADE can give new results. 

4. The presented asymptotic procedure (based on the equation for vorticity) is re- 
markably simple. However it has one disadvantage: modified pressure in ADE (|3.24|) and 
(|3.25p remains undefined (in terms of physical parameters involved). Therefore the im- 
portant next step is to repeat the same asymptotic procedure for the Euler equations, 
from the very beginning. The expressions for modified pressure will allow one to calculate 
the forces exerted on walls, particles, etc. 

5. The notion of a drift is actively used in this paper. It is known that a drift velocity 
can appear from Lagrangian, Eulerian, or hybrid (Euler-Lagrange) considerations; look 



for the discussion related to these options in Vladimirov (2012) 



6. The relation between high-frequency asymptotic solutions and small-amplitude 
asymptotic solutions were studied in Vladimirov (2011)] in particular, for purely periodic 
non-modulated solutions these classes can be isomorphic to each other (any solution from 
one class can be transformed into a solution from another class). 

7. The mathematical justification of the derived averaged equation by the estima- 
tion of an error in the original governing equation can be performed in a similar to 
Vladimirov (2010)| |Vladimirov (2011)| way. 



8. One can also derive the higher approximations of ADE as it has been done in 
Vladimirov (2010)[ |Vladimirov (2011)| They can be especially useful for the study of 



motions with Vq = 0, see Vladimirov (2010) 



9. Viscosity and diffusivity can be incorporated into our asymptotic theory However, 
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viscous and diffusion terms produce additional to ()3.3|) independent small parameters 
and the result will essentially depend on their order in terms of e. In particular, the scale 
of slow time s (|3.5p can be different. 

10. ADE is valid not only for infinite flows but also for finite flow domains, including 
the boundaries that radiates sound. 

The author is grateful to Profs. A.D.D.Craik, S.Leibovich, H.K.Moffatt, and N.Riley 
for helpful discussions. 
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